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ABSTRACT 

A  brief  review  of  the  operation  of  vibratory  gyroscopes 
and  the  status  of  commercial  micromachined  gyroscopes  is 
presented  in  the  introduction.  Following  the  introduction, 
an  alternative  design  approach  employing  multiple, 
coupled  gyroscopes  per  angular  axis  is  examined  as  a 
method  to  improve  performance  and  redundancy  in 
angular  rate  sensors.  The  results  arc  from  simulations  and 
explore  the  effects  of  non-identical  gyroscopes,  variations 
in  the  driving  frequency,  and  coupling  on  synchronization 
in  the  array.  A  novel  operating  approach  is  presented  that 
requires  the  array  of  gyroscopes  to  be  synchronized  in 
phase  and  frequency  to  improve  the  detection  of  sense 
axis  displacement  and  to  utilize  one  amplitude  demodulator 
for  the  entire  array  of  gyroscopes. 

1.  INTRODUCTION 

Vibratory  gyroscopes  are  designed  to  determine  the 
angular  rale  of  a  rotating  object.  Typically,  vibrating 
gyroscopes  can  be  modeled  by  the  mass  spring  system 
shown  below  in  Figure  1 . 


Figure  1.  Model  for  vibratory  gyroscope 

A  mass  spring  system  used  to  model  difTcreni  types  of  vibratory 

gyroscopes. 

The  operating  principle  for  all  vibrating  gyroscopes  is  the 
transfer  of  energy  from  one  vibrating  mode  to  another  by 
the  Coriolis  force.  The  Coriolis  force  is  an  inertial  force 
exerted  on  a  moving  body  in  a  rotating  reference  frame. 
The  gyroscope  in  Figure  1  has  two  orthogonal  excitation 
modes  with  restoring  coefficients  ky  and  k*.  Damping 
coefficents  along  each  axis  (C^  and  C,)  arc  included  in  the 


model.  The  y-axis  mode  (drive)  is  excited  by  a  reference 
drive  force  and  the  z^axis  mode  (sense)  is  excited  by  the 
Coriolis  force  (F„ ).  Displacements  along  the  sense  or  z 
axis  are  proportional  to  the  applied  angular  rate  O,,  if  the 
motion  of  the  drive  axis  is  known.  The  following  equations 


model  this  system. 

~  ~^yy  “  y  dHve  QXIS 

(la) 

mz  =  -  C,i  +  sense  axis 

(lb) 

=  Coriolis  Force 

(Ic) 

=|  2m  X  V  I 

The  critical  design  i.s.sues  for  micromachined  vibratory 
gyroscopes  are  detecting  sense  axi.s  displacements  that  are 

orders  of  magnitude  less  than  the  drive  amplitude  and 
keeping  the  drive  and  sense  motions  uncoupled.  The 
excitation  of  the  drive  mode  and  displacements  along  the 
sense  axis  arc  typically  measured  using  electrostatic 
sensors  by  applying  a  voltage  or  measuring  changes  in 
capacitance.  The  amplitude  of  the  drive  axis  motion  is 
typically  limited  to  the  1-50  micron  range  depending  on  the 
etching  and  manufacturing  process  used.  Typical  sense 
axis  displacements  arc  orders  of  magnitude  less  than  drive 
axis  displacements  and  require  sensitive  transducers  to 
measure  the  small  displacements.  This  paper  will  focus  on 
low  cost  surface  micromachined  gyroscopes  with 
capacitive  sense  axis  detection.  While  capacitive  sensors 
are  not  ideal  for  micromachined  applications,  they  can  be 
implemented  using  standard  silicon  processkig 
techniques.  Micromachined  gyroscopes  with  bias  drifts  in 
the  10-50®/h  have  been  produced  on  a  prototype  basis  [1]. 
An  application  for  low  cost  gyroscopes  with  a  10®/h  bias 
drift  is  in  the  automotive  industry  for  roll  sensors. 
Gyroscopes  with  a  bias  drift  less  0.l°/h  are  suitable  for 
inertial  navigation  applications. 

2,  SIMULATION  MODEL 

The  equations  in  (1)  can  be  put  into  a  more  general  form  if 
an  amplitude  modulated  sinusoidal  driving  force  is  used  to 
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represent  the  Coriolis  term  in  the  sense  axis  equation  (lb). 
This  simplification  requires  the  drive  axis  displacement 
follow  a  harmonic  motion  of  the  form  cosCw^^t).  This  can 
be  accomplished  by  locking  the  drive  axis  displacement  to 
a  harmonic  reference  signal  in  a  gain  control  loop. 
Equation  (2)  Ascribes  the  sense  axis  motion  (z*axis  in 
figure  1)  of  the  gyroscope  using  a  sinusoidally  driven, 
linearly  damped,  harmonic  oscillator.  Equation  (3)  casts 
equation  (2)  in  more  general  terms. 


mz'VCi^rkz^  A:n(/)  Xsin  (w^/)  (2) 

z  +  2yz  +  il(t)Xsin  (w^/) 


In  equation  (3)  y  is  C/2m,  X  is  a  constant  that  includes  the 
drive  amplitude,  is  the  harmonic  drive  frequency,  w,  is 
the  resonant  frequency  for  the  sense  axis  mode,  and  ft(t) 
is  the  angular  rotation  rate  about  the  x-axis.  The  general 
solution  for  z(t)  can  be  found  and  it  describes  the 
displacement  of  the  sense  axis  of  the  gyroscope  [2]. 
Equations  (4)  and  (5)  follow  from  the  general  solution  in 
the  underdamped  case  where  y  <  describe  the 

maximum  amplitude  of  the  sense  vibration  mode  and  (p  is 
the  phase  lag  of  the  gyroscope  with  respect  to  the  driving 
frequency  and  the  gyroscope  sense  axis  resonant 
frequency  (w^). 


An(/) 


1  M'r  K<2 


0  =  arctan 


w. 


Vt' 


,  0  ^0  < ;r 


(4) 

(5) 


These  equations  describe  the  necessary  conditions  for  the 
gyroscopes  in  the  array  to  synchronize  and  will  be  used 
extensive  to  qualitatively  describe  the  results  of  the 
simulations. 


Figure  2b.  Plot  of  Phase  Lag,  ^  -  This  figure  shows  that  for 
Wj  >  wj  the  gyroscope  will  follow  the  drive  frequency  with 
roughly  zero  phase  lag.  For  the  sense  mode 

follows  the  drive  frequency  by  For  >  w,  the  sense 
mode  follows  by  K.  The  region  of  most  phase  mismatch 
will  be  near  w,  and  will  be  worse  for  high  Q. 

The  main  result  is  that  for  driving  frequencies  less  that  w, 
the  sense  mode  oscillates  approximately  in  phase  with  the 
drive  frequency.  The  design  of  an  array  of  non-identical 
gyroscopes  will  involve  choosing  a  Q  for  the  sense  axis 
that  minimizes  pha.se  mismatch  in  the  array  and  maximizes 
the  amplitude  response  to  the  Coriolis  driving  force.  The 
simulation  results  presented  in  the  next  section  explore 
some  of  the  tradeoffs  in  such  a  design. 

3.  SIIVfUT ATION  RESULTS 


Q  is  equal  to  wj(2y)  and  is  called  the  “quality  factor”  of 
the  sense  mode  of  the  gyroscope.  Plots  of  equations  (4) 
and  (5)  arc  shown  in  figure  2. 


Mix  Amplitude  ai  hinctlon  Q 


Figure  2a.  Plot  of  Zmax  -  This  figure  shows  ihc  sharpening  of 
the  amplitude  response  as  approaches  w,  and  Q  increases. 


For  all  of  the  simulations  an  array  of  one  hundred 
gyroscopes  in  a  ten  by  ten  two  dimensional  array  was 
used.  The  mean  values  for  the  mass  (m).  linear  restoring 
forcc(kl),  damping  (C).  drive  resonance  (w,),  applied 
angular  rate  (fl).  drive  amplitude  (Aj) ,  and  drive  frequency 
(wj)  are  shown  in  table  1 . 


mass 

lc-9 

kg  

Kl 

1 

N/mctcr 

C 

3.1623C-7 

N  t/mcicr 

31623 

rad/sec 

a 

0.175 

rad/sec 

Ai 

lOc.6 

meters 

JHa _ 

31623 

rad/sec 

Table  I.  Simulation  Parameters  with  no  variation 
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The  simulations  begin  with  an  “ideal’*  case  where  the 
gyroscopes  are  identical  and  each  gyroscope  is  driven 
exactly  at  After  the  initial  plots,  these  constraints  are 
removed  and  the  mass  is  allowed  to  vary  by  le-9±2.5e-ll 
kg  with  equal  probability.  This  creates  a  variation  in  the 
resonant  frequency  of  the  drive  axis  of  5033Hz  ±  62.5  Hz. 
For  many  of  the  simulations  wc  will  look  at  drive 
frequencies  equal  to  the  resonant  frequency  of  the  sense 
axis  and  at  4966  and  5100  Hz.  There  drive  frequencies  lie 
outside  the  maximum  variation  in  the  sense  frequency 
introduced  by  varying  the  mass  of  the  gyroscope  by  2.5%. 
For  all  the  simulations  a  random  set  of  initial  conditions 
was  imposed  on  the  system.  The  initial  displacements  of 
the  sense  mode  for  each  gyroscope  was  randomly 
distributed  between  0  ±  l.5e-8  meters  with  equal 
probability.  The  value  of  l.5e-8  is  roughly  the  maximum 
displacement  for  a  single  gyroscope  with  the  parameters 
given  above.  The  sense  mode  was  also  given  a  random 
velocity  between  0  ±  I.Oe-4  meters/sec  with  equal 
probability.  Initial  velocities  larger  than  this  tended  to 
destabilize  the  system.  These  random  initial  conditions 
place  the  array  in  very  disordered  state.  These  simulations 
do  not  reflect  the  expected  performance  of  the  gyroscopes. 
They  were  completed  to  explore  and  optimize  the 
synchronization  performance  of  the  array  of  gyroscopes. 
A  constant  external  angular  rate  of  0.175  rad/sec  was 
applied  in  all  simulations  and  the  time  for  the  .summed 
displacements  of  each  gyro  to  reach  a  constant  amplitude 
was  used  to  measure  synchronization. 

Figures  3a -c  illustrate  the  initial  conditions  imposed  on  the 
array  and  the  summed  output  and  amplitude  demodulated 
signal  from  the  array  of  one  hundred  gyroscopes. 


Figure  3a-c.  Time  scrict  and  Demodulation.  All  figures  use  an 
ideal  model  with  identical  gyroscopes  and  drive  frequency  at 
resonance  (5033Hz).  3(a)  shows  the  effect  of  initial  random 

conditions  from  t*0  to  t-lOO’Ts,  where  Ts  *  20*5c-7  seconds. 
3(b)  shows  the  cxcelicnl  phase  and  frequency  synchronization  from 
t“1900*Ts  to  2000*Ts.  3(c)  shows  the  increasing  amplitude  of  the 
summed  output  from  each  gyroscope  and  the  demodulated  signal. 
This  demodulation  method  is  only  possible  for  phase  and  frequency 
synchronized  gyroscopes. 

Figure  4  shows  the  results  from  an  “ideal”  system  with  identical 
gyroscopes  and  no  variation  in  the  drive  frequency.  A  coupling 
term  is  also  introduced  in  the  model  which  has  the  form 
where  Z  is  the  sense  mode  displacement 
and  kappa  is  the  coupling  strength.  A  kappa  of  zero 
corresponds  to  an  uncoupled  system.  In  the  simulations,  a 
nearest  neighbor  coupling  method  was  employed.  For 
gyroscopes  in  the  center  of  the  array  the  coupling  term  includes 
its  four  neighbors  and  the  coupling  for  edge  elements  contains 
fewer  terms. 

The  coupling  term  could  be  implemented  in  gyroscopes  through  a 
variable  resistor  connecting  the  sense  axis  capacitors  of  nearest 
neighbors  in  the  array.  The  coupling  term  was  modeling  in  these 
simulations  to  see  if  it  improved  the  amplitude  response  or 
decreased  the  synchronization  time  for  the  summed  gyro  output 
signal.  Figure  4  shows  the  effect  of  coupling  on  the  ideal  system 
with  identical  gyroscopes.  No  improvement  in  the  amplitude 
response  or  synchronization  time  is  achieved  using  coupling. 


Figure  4.  Identical  Gyroscopes  driven  at  resonance  (5033 
Hz).  This  plot  shows  that  coupling  decreases  the  amplitude 
response  for  an  ideal  array  of  gyroscopes. 
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in  Hgurc  5  a  more  realistic  system  that  includes  a  random  2.5% 
variation  in  the  mass  ts  examined.  This  mass  variation  changes 
the  resonant  frequency  of  the  sense  mode  between  5033  ±  62.5 
Hz.  This  system  now  docs  not  have  a  single  resonant  frequency 
to  exploit.  Three  values  for  the  drive  frequency  were  chosen  to 
examine  the  response  of  the  array.  They  are  at  resonance,  and 
slightly  above  and  below  the  frequencies  variation  possible  for 
the  sense  axis  mode. 


Figure  5.  Non-Identieal  Gyroseopea  and  Effect  of  Drive 
Frequency. 


The  effect  of  coupling  does  not  change  the  general  response  of 
the  array  of  gyroscopes  but  it  does  decrease  the  maximum 
amplitude  of  the  summed  gyroscope  output  signal.  The  final 
case  that  will  be  considered  is  minor  variations  in  the  drive 
frequency  Wj. 


Figure  7.  Non-Identical  Cyroi  with  Variation  in  Drive 
Frequency 


The  system  driven  at  resonance  now  exhibits  the  worst 
amplitude  and  synchronization  petfoimance.  This  is  due  to  the 
maximum  slope  change  in  the  phase  lag  lesponse  near  the  mean 
resonance  frequency  (see  figure  2).  The  best  amplitude  response 
is  for  a  driving  frequency  above  resonance  (5100  Hz).  This  can 
be  explained  qualitatively  by  realizing  that  in  this  regime  all  the 
gyroscopes  meet  the  condition  w,  <  and  thus  have 
approximately  a  ti  phase  lag  with  respect  to  the  drive  frequency. 
This  phase  lag  does  not  affect  the  demodulation  of  the  summed 
output  from  the  array  since  the  gyroscopes  arc  synchronized  to 
the  drive  frequency  with  a  constant  phase  lag. 

In  figure  6  a  similar  plot  to  figure  5  is  done  with  the  coupling 
strength  kappa  fixed  at  O.I.  This  figure  is  very  similar  to  the 
uncoupled  .system. 


Figure  6.  Non-Identical  Gyroicopes  and  Effect  of  Drive 
Frequency 


The  top  three  lines  in  figure  7  use  a  drive  frequency  equal  to  the 
mean  resonance  frequency  of  the  array  (5033Hz)  and  a  variation 
in  the  drive  axis  of  only  0.5  rad/scc.  The  lower  three  figures  use  a 
drive  frequency  greater  than  the  mean  sense  resonant  frequency 
and  the  variation  in  the  drive  is  12.5  rad'scc.  In  both  cases  the 
coupled  system  performs  better  than  the  uncoupled  system. 

4.  SUMMARY 

This  paper  demonstrates  that  for  non-identical  arrays  of 
gyroscopes  with  variations  in  the  drive  frequency  coupling  can 
increase  the  synchronization  in  the  array.  The  simulations  and 
solutions  for  a  single  damped,  driven  harmonic  oscillator  also 
indicate  that  a  driving  frequency  above  the  mean  variaiion  in 
resonant  mode  frequencies  produces  the  greatest  synchronization 
in  the  array.  Future  work  will  explore  the  effect  of  noise  and  a 
lime  varying  angular  velocity.  Ti  will  also  he  necessary  to 
characterize  the  response  of  the  summed  array  displacement 
output  to  the  applied  angular  velocity. 
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